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ABSTRACT 

The chemistry in the central regions of galaxies is expected to be strongly 
influenced by their nuclear activity. To find the best tracers of nuclear activity is 
of key importance to understand the processes taking place in the most obscured 
regions of galactic nuclei. In this work we present multi-line observations of CS, 
C^^S, HNCO and C^^O in a sample of 11 bright galaxies prototypical for different 
types of activity. The ^^S/^^S isotopic ratio is ~ 10, supporting the idea of an 
isotopical ^^S enrichment due to massive star formation in the nuclear regions of 
galaxies. Although C^^S and C^^'S do not seem to be significantly affected by the 
activity type, the HNCO abundance appears highly contrasted among starburst. 
We observed HNCO abundance variations of nearly two orders of magnitude. 
The HNCO molecule is shown to be a good tracer of the amount of molecular 
material fueling the starburst and therefore can be used as a diagnostics of the 
evolutionary state of a nuclear starburst. 



Subject headings: ISM: molecules — galaxies: abundances — galaxies: ISM — 
galaxies: starburst — galaxies: Seyfert — radio lines: galaxies 
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Introduction 



The nuclear regions of active galaxies can harbor different energetic phenomena such as 
starbursts (SBs) and Active Galactic Nuclei (AGN) which are responsible for producing the 
bright emission stemming from their central regions. The heating of the interstellar medium 
(ISM) in nuclei with different type of activity is expected to be dominated by a variety of 
mechanisms, such as UV, X-rays, and/ or shocks, which will ultimately drive their observed 
chemical richness ( iMartfn et al.ll2006bl . and references therein). Disentangling which of this 
mechanisms is the main power source of galactic nuclei is complicated by the high obscuration 
affecting these objects. This fact becomes critical in the case of Ultra Luminous Galaxies 
(ULIRGs) and high-z sources. Therefore the observation of the dense molecular material has 
become an essential tool to get an insight in the evolution and classification of the nuclear 
activity in galaxies. 

Finding appropriate molecular species to accurately trace each of these heating mech- 
anisms is crucial to stablish the nature of the central engine. Molecular species like HCN 
and HCO^, and in particular the ratios HCN/CO and HCN/HCO^, have been claimed to 



Kohno et al. 


1999: 


Kohno 


2005; 


Krips et al. 


2008) 



are though t to show enhanced abundances in extremely irradiated environments (lAalto et al. 



2OO21 . 120071 ). A dichotomy still remains concerning the different evolutionary state of nuclear 
starburst. While at an early stage, the heating of the ISM is thought to be dominated by 
shocks affecting the molecular clouds fueling the starburst, the late stages of starburst are 
vastly dominated by the UV radiation originated in the newly formed massive stars. This 
scenario has been inferred from the extensive observation of the prototy pe sources NGC 253 



and M 82 in a number o f key molecular tracers s uch as SO2, II2S, OCS (IMartin et al.l 12003 



2005h . HOC+, and C0+ dPuente et aPboodbood l Unfortunately, these tracers are generally 
too faint to be detected, but for the brightest prototypical nearby galaxies. 

In a recent observational study carried out towards a sample of mo lecular clouds dom i- 
nated by different heating mechanisms within the Galactic center region (IMartm et al.ll2008l ). 
we have found that the CS/HNCO abundance ratio is highly contrasted between molecular 
clouds illuminated by the UV radiation from massive star clusters and the giant molecular 
cloud complexes shielded from photodissociation and mostly heated by shocks. The origin 
of the large changes in the CS/HNCO ratio is due to the enhancement of CS in phot on- 
dominated regions (PDRs) through reactions involving S"*" (jSternberg fc Dalgarndll995f ) as 
opposed to the highly photodissociable HNCO, which is destroyed by the UV photons. Both 
observations of photodissociation region s within our Galaxy and photochemical models show 
the suitability of CS as a PDR tracer (IGoicoechea et al.ll2006l ). In addition, HNCO is en- 
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hanced in the presence of shocks due to its injection in the gas phase from the grain mantles 
( iZinchenko et aLl 120001 ) . This paper presents a follow up study of the CS/HNCO ratio to- 
wards a sample of prototype galaxies wit h different nuclear activity. Although CS has been 
extensively observed in external ga laxies (Baan et al. 20081), HNCO had only bee n detected 
in four galaxies prior to this work (jNguyen-Q-Rieu et al.lll99ll : IWang et al.ll2004j ). 



2. Observations and Results 

Observations were carried out with the IRAM 30 m telescope on Pico Veleta, Spain, 
during three observing periods from summer 2005 through summer 2007. Table [T] shows the 
sample of galaxies with their coordinates and the rest frequencies of the observed molecular 
transitions. Observations were performed in symmetrical wobbler switched mode with a 
frequency of 0.5 Hz and a beam throw of 4' in azimuth. As spectrometers we used the 
512 X 1 MHz filter banks for the 3 mm transitions and the 256 x 4 MHz filter banks for those 
at 2 and 1.3 mm. Pointing accuracy was estimated to be of ~ 3" from frequent continuum 
cross scans on nearby pointing sources. Data were calibrated using the standard dual load 
system and main beam temperatures were obtained as Tub = {Fcs/Bcs)Tl, where i?efr is 
tabulated in Tableland F^s are 0.95, 0.93, and 0.91 for 3, 2, and 1.2mm, respectively. 
Calibration uncertainties are estimated to be < 15%. 

The observed spectra for each source are shown in Fig. [H The HNCO 5 — 4 and 10 — 9 
lines were observed simultaneously with the C^^O 1-0 and 2 — 1 lines, respectively. This 
allows a reliable determination of the HNCO relative abundance, as no uncertainties due 
to pointing drifts and/or calibration are involved. For some of the sources with the larger 
observed linewidths, the HNCO and C^^O lines appear moderately blended. Such is the case 
of NGC7469, NGC1068, and Arp220. Gaussian shapes have been fitted to the observed 
profiles. Table [2] compiles the derived fitted parameters. For the non-detected transitions, 
upper limits are tabulated where it was considered a 3a level and an approximate linewidth 
similar to those of the detected CS and HNCO transitions, or that of C^^O in the case of 
NGC5194. 

Rotational temperatures and molecular column densities have been derived assuming 
local thermodynamic equilibrium (LTE) and optically thin emission. As derived for CS 
J = 3 — 2 in NGC253 ( Martin et al.ll2005l ) we only might expect the emission of CS J = 5 — 4 
to be significantly optically thick. However, our observation of C34S will be less affected 
by the opacity effects. We have accounted for the dilution of the source in the telescope 
beam by assuming an averaged emission extent of 10" for all sources. This simplification 
might affect the determination of the column densities for sources much smaller than 10" 
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(see iMartm et al.ll2006bl ). but the abundance ratios between species dealt with in this paper 
(Sect. will not be significantly affected. Thus, the main source of uncertainty will be due 
to possible different emission extents between molecules. Table [3] shows the column density 
and rotational temperature derived for every source and species. From the LTE analysis, 
we derive C^^O rotational temperatures (Trot) of 3 — 7K for most of the galaxies except for 
M 82 where we determine Ti.ot=15 K. The temperatures derived from the HNCO line emission 
of 4 — 18 K are a factor of ~ 2 higher than those derived from C^^O. Whenever only one 
transition was observed, the average rotational temperature of 10 K was assumed. 

We have also included th e available observations towards NGC 253 to derive its physica l 
parameters shown in Table [3] ( Harrison et al.lll999l : iMartm et al.ll2005l : iMartm et al.ll2006bl ). 
Note that the C^^O transitions in NGC 253 were observed ^ 10^^ north of the position 
shown in Tabled! For NGC 4945 we used the observations by IWang et al.l (120041 ). Not all 
sources in the sample were observed in th e CS J = 5 — 4 transitioii . The upper limit to the 



emission to this transition in NGC 6946 (IMauersberger et al.l 1 19891 ) allows to constrain the 
CS column density to < 6.2 x lO^^cm"^. We have thus used the CS J = 3 — 2 transitions 



available on NGC 6946 and NGC 5194 (IMauersberger et al.l Il989l ) to estimate the column 
densities presented in Table [ SI Note that for NGC 6946 we used the observed offset (Aa = 



10", A(5 = —10") position in IMauersberger et al.l (119891 ). which is located just ~ 1" away 
from our nominal position. Similarly, their NGC 5194 (0", 0") position is ^ 4" away from 
our observed position. 

The critical densities for CS and HNCO transitions observed in this works are expected 
to be very simi lar since the Einste in coefficients and coUisional cross sections are estimated 
to be similar (ISchoier et al.l l2005l ). Thus, the HNCO/CS ratio should refiect a reliable 
abundances ratio. In any case, the possible uncertainties introduced by shghtly different 
critical densities will never account for the large difference in the abundance ratios measured 
in different galaxies. The differences in the derived molecular abundances might also be 
affected by molecular excitation differences of each individual source. However the detailed 
multitransition study and LVG model i ng of dense gas tracers such as HCN and HCO+ 
in a sample of galaxies by iKrips et al.l (120081 ) . shows how for several molecular transitions 
with similar critical densities, the derived abundances from the line intensities are really 
determined by real differences in abundances and not by excitation effects. The low excitation 
temperatures derived from C^^O indicates that a substantial fraction of the column density 
will be dominated by relatively low density gas. Since in our diagnostic diagrams in Fig. [3] 
are normalized to the total column densities derived from C^^O, changes in the H^ total 
column densities will not change the main conclusions of the paper. 
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3. Selected sources 

The sample of sources observed in this paper has been selected among the brightest 
prototypes for the different types of nuclear activity. Their distances and corresponding 
linear scales are shown Table [TJ The sources can be grouped as follows. 



3.1. Starburst Galaxies (SBG) 



M82 and NGC 253 are the strongest and richest extragalactic molecular sources. These 
are the archetypes of starburst galaxies housing the two brightest extragalactic IRAS sources 
(ISoifer et al.lll989l ). M82, as opposed to the young starburst in NGC 253, is claimed to be 
the prototype of evolved starburst. This is supported by its ob served chemistry characterized 



by low abundance of complex molecules such as SiO, CH3OH (IMauersberger &: Henkellll993 



Martin et al.l l2006bl ) and large abundance of r nolecular ions like CO"*", HOC"*", and H3O 



which are expected to be enhanced in PDRs (IFuente et al.l l2005l . l2006l : Ivan der Tak et al. 



20081 ). In this work we did not observe the central position towards M82, but a position in 
the north-east molecular complex (hereafter M82*) where most of the photodissociation re 



gions are located as observed in the HCO emission high-resolution maps (iGarcia-Burillo et al. 



2OO2I ). Majfei 2, is another well studied starburst spiral galaxy with a high molecular gas con- 
centration i n its nuclear region and showing traces of a tidal interaction with a dwarf compan- 
ion galaxy (iHurt et al.lll996l : lMason fc Wilsonll2004j ) . M 55 is a nearby fa ce on barred galaxy 
unde r going strong nuclea r starburst likely due to gas inflow along its bar (jPetitpas &: Wilson 



1998 



Talbot et al.l Il979l ) . NGC 694 6 s hows a moderate starburst in its nucleus a lso claimed 



to be related to the presence of a bar (ITurner fc Holll983l : ISchinnerer et al.ll2007l ). 



3.2. Active Galactic Nuclei (AGN) 



NGC 1068 is a nearby luminous infrared galaxy (LIRG) prototype of a Seyfert 2 nucleus 
jSf) . The AGN is enclosed by a circumn uclear starburst ring with 14" 



( lAntonucci &: Mi 



lei 



198^ 



(~ 1 kpc) radius (IMyers fc Scovillelll987l : ISchinnerer et al.ll2000l ). Thus, we have observed two 
positions in NGC 1068. One towards the central AGN and an offset position 0", —14" towards 
a peak of emission within the circumnuclear ring, hereafter NGC 1068*, mostly tracing the 
starburst. NGC 4945, being one of the thre e brightest I RAS point sources, has been the 
target of some a detailed molecular study by IWang et al.l (120041). This nearly edge-on spira l 
galaxy harbor a heavily obscured Seyf ert 2 nucleus (iBraatz et al.lll997l : iMaiolino et al.lll999l ) 
also surrounded by a starburst ring (IMarconi et al.ll2000l ). NGC 5194 (M51) is a nearby 
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grand-design spiral galaxy with a Seyfert 2 nucleus. The proximity of this galaxy as well as 
its nearl y face-on orientation has made it the target of numerous large scale multi- wavelength 
studies ( IScoville &: Young Il983l : IScoville et al.l I2OOII : ISchuster et al.l 120071 ) . Its interaction 
with the companion g alaxy NGC5195 seena s to have been the orig in of a past period of 
intense star formation (IThronson et al.lll99ll : iGreenawalt et al.lll998l ). NGC 7469, the only 
example of Sey fert 1 nucleus in ou r sample, also shows a circumnuclear star- forming ring of 
3" in diameter (j Wilson et al.lll99ll ). 



3.3. UltraLuminous Infrared Galaxies (ULIRG) 



Arp 220 is one of the best studied nearby ULIRGs (IDopita et al.ll2005r) w ith more than a 
95% of its luminosity radiated at IR/submm wavelengths (ISanders et al.ll2003l ). At a distance 
of D~ 77Mpc {z = 0.018) and star formation rates of 3OOM0yr~^, Arp 220 is often referred 
as a nearby template of the luminous star-forming galaxies at high-z. Arp 220 is the result 
of an advanced merger syst em with two nuclei separated by ~ 1" (~ 3 00 pc) as seen from 
near-IR and radio images (IBaan fc Haschickl Il995l : IScoville et al.l 120001) with an enormou s 
concentration of molecular gas within its central region (a few 10^ Mg. lSakamoto et al.lll999l ). 



3.4. "Normal" 



We refer to IC 342 is as a "normal" galaxy as being the closest spiral galaxies resem- 
bling our own Galaxy. It shows a mini spiral structure meeting at the inner rnolecular ring 
surro unding the central stellar cluster (jlshizuki et al.l Il990l : iBoker et al.l Il997l : iHelfer et al. 
2OO3I ) . Chemical differences have been observed between the gas in the nuclear region, af- 
fected by the intense star formation through the radiation from the c entral cluster, and the 



molecular raateria l in the spiral structure, mainly affected by shocks (iMeier &: Turnerl 12005 



Usero et al.ll2006l ). In this work we have observed two positions, one at the central region 
and the other at the offset position (+5", +15"), hereafter IC 342* , located just in one of the 
arms of the minispiral. 
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Discussion 



4.1. C^^^S and C^^S 



We have derived the abundances of CS, C S, and HNCO relative to H2. In order 
to estimate the H2 column density we have used the obser ved C^^O emission, assuming 



a 0/ abundance ratio of 150, as derived in NGC253 flHarrison et al.l Il999l ) and the 
standard CO/H2 ratio of 10~^. 

The CS fractional abundance shows a narrow range of values within 1 — 3 x 10^^. 
Similarly the C^^S abundances are concentrated in the range ofl — 5x10^^*^. There are only 
three exceptions to this trend, namely IC342*, Arp220, and NGC5194. The abundances 
of CS observed towards IC342* is a factor of 6 below the average in the sample, while the 
C^^S abundance is close to the average value. This is explaine d by the highly concentrated 



CS emission towards the nuclear region (IMeier fc Turner! |2005| ) and the almost twice smaller 
beam size at the CS 5 — 4 transition frequency compared to that at the C^^S 3 — 2 line (see 
Tabled]). We did not detect neither CS nor C^^S emission towards NGC5194. While the 
C^^S upper limit is still within the average value obtained for the other galaxies, the limit 
we derive in the main isotopologue is a factor of 3 below the sample average. On the other 
hand, the galaxy Arp 220 is observed to have an abundance a factor of ~ 5 over the average 
in both isotopologues. This points out a real overabundance of CS towards this ULIRG. 

Fig. [2] shows the C^^S/C'^^S abundance ratio versus the fractional abundance of CS. 
This ratio is calculated with the 1 mm CS 5 — 4 and the 2 mm C^*^ 3 — 2 transitions which 
makes it dependent on the assumed rotational temperature to derive the column densities, 
and to a lesser extent (up to a 30% for the smallest sources) on the considered source size. 
However, the ratio derived for NGC4945 and NGC 253 agree wi t hin a facto r of 1.5 with 



the va lues derived from the complete study of CS by IWang et al.l ( 120041 ) and iMartm et al. 



( l2005l ). respectively. Although the effects of line opacity have not been taken into account, we 
can consider the derived C^^S/C^^'S abundance ratio as a rough representation of the ^2/34g 
isotopic ratio. It appears that almost all the sources in the sample have ratios around ~ 10, 
in agreement with those previously derived for NGC4945 and NGC 253, and well below the 
value of ~ 24 measured in the Galactic disk (IChin et al.lll996l ). This result supports the idea 
of an ^^S overproduction in the nuclei of galaxies by massive stars. Furthermore, we observe 
this ratio to be independent of the type of nuclear activity, both in AGN and SB dominated 
galaxies. 
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4.2. CS versus HNCO abundances 

In Fig. [3l we represent the fractional abundances of CS and C^^S versus those of HNCO. 
We observe that most of the galaxies in the sample are observed with HNCO fractional abun- 
dances ranging from 10~^ to 10~®, but for two sources, namely NGC5194 and M82, showing 
remarkable low HNCO abundances. For these sources we derived fractional abundance below 
the average derived from the rest of galaxies in our sample by a factor of ~ 3 and ~ 20, 
respectively. 

We do not observe major differences in the derived abundances among the AGNs and 
SBGs in our sample. This homogeneity suggests that the type of nuclear activity does not 
play a major role at large scales in the production/destruction of these species in the bulk 
of the molecular gas. This effect is illustrated by the two positions observed in NGC 1068, 
where the HNCO and ^^CS abundances decrease by less than a factor of two from the nuclear 
AGN to the SB ring. Of course we have to bear in mind that the two regions are not fully 
resolved by the beam of the IRAM 30 m telescope at these frequencies. 



If, as suggested by the observations in the Galactic Center region (IMartm et al.ll2008l ). 
the abundance of these species is mostly driven by the presence of shocks and strong UV 
radiation fields in the ISM, then the observed differences are likely linked to star formation 
activity in these galaxies. Significant differences are observed within the galaxies IC 342 and 
Maffei2. In IC342, the HNCO/C^^O ratio increases towards the position in the minispiral 
and decreases in the central region where ISM is illuminated by the UV radiation from the 
central clusters. These chemical differences, smoothed to the 10" — 20" singl e-dish beam size, 



are m uch more contrasted in the high-resolution interferometric maps by (IMeier fc Turner 



20051 ) . Although not spatially, we resolve the nuclear region of Maffei2 by the tw o velocity 



components. As seen from the high-resolution maps of ^^CO (IMeier et al.ll2008l ) the lower 
velocity component ( — lOOkms"^) is associated to the nuclear and northeastern molecular 
arm, while the higher velocity component (Okms~^) will be mostly associated to the nuclear 
and southwest molecular arm. We clearly see a difference in the chemistry of both compo- 
nents with an increase of the HNCO abundance towards the southwest molecular arm of a 
factor of ~ 2 as well as a decrease of the CS abundance of a factor of 3 — 4. Overall, the 
HNCO/CS ratio increases by almost an order of magnitude. 

We can compare the ratios r neasured in galaxie s with the abundances derived in the 



central region of the Milky Way (IMartm et al.ll2008l ). The vast majority of galaxies show 



HNCO abundances intermediate between those found in the shock dominated giant molec- 
ular clouds (GMCs) and the Galactic PDR prototypes (see dashed lines in Fig [3]). Two 
sources in our sample, namely M83 and IC342* clearly fall in the range of abundances 
found in GMCs. On the contrary, M82 represents the other extreme where HNCO emis- 
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sion is utterly vanished. The upper hmit to the detection derived after a deep 5.5 hours 
integration on the HNCO 5 — 4 is an o rder of magnitude lower th an the previous limit 
to the detection in the central position (iNguyen-Q-Rieu et al.lll99ll ). but still HNCO re- 
mains undetected. Interferomet ric maps with the BIM A array shows no emission anywhere 
in the nuclear region of M82 (ITurner fc Meierl 120081 ). This non-detection of HNCO sup- 
ports the idea that M 82 is highly pervaded by a strong UV photodissociating radiation and 
depleted of dense molecular a s also suggested by the low abundance o f CH3OH and NH3 
( iMauersberger fc Henkell Il993l : iTakano et al.l l2002l : iMartm et al.l l2006al ) . We would expect 
to observe a hig her abundance of CS in this NE position where photodissociation is claimed 
to be dominant (iGarcia-Burillo et al.ll2002l ). However, we observe just a slightly higher CS 
abundance than in the whole sample, and not even that if we consider the isotopologue 
C^^S. Altogether with M82, the upper limit to the HNCO abundance derived for NGC5194 
are similar to those derived in the Galactic Center PDRs. Most of the molecular studies 



of NGC5194 so far are restricted to the emission of CO and HCN (IKohno et al.lll996l . and 
references therein). The results from this work point out the possible similarities between 
the chemistries of NGC5194 and M82, both with very low abundances of HNCO as com- 
pared with the rest of galaxies in this sample. This result agrees with the idea of NGC5194 
being in a post starburst stage after th e star formation event triggered by its interaction 
with NGC5195 JCreenawalt et al.lll998 [). Additionally, the HNCO and CS abundances in 
NGC253 and NGC1068*, a factor of 3 below the average, as well as the significantly high 
C^^S abundance, could be explained by photodissociation. This would suggest that the UV 
radiation from the evolved SB plays a significant role in the chemistry of the nuclear ISM in 
these galaxies. Therefore NGC253 would contain a significant PDR molecular component, 
similar to what is observed in M82. However, the amounts of dense molecular gas shielded 
from the UV radiation is still one order of magnitude larger than that observed in M82. 

Our observations of HNCO confirm that this molecule is one of the best tracers to study 
the evolutionary state of nuclear starbursts in galaxies because of the high contrast shown 
between the UV pervaded evolved starbursts and the dense and well shielded molecular ma- 
terial in the early stages of star formation. This molecule shows a hi gher abundance variation 
in galaxies than oth er dense molecular gas tracer such as methanol (iHiittemeister et al.lll997l : 
Martin et al.ll2006al ). On the other hand, CS does not show a significant variation among the 
observed sources. This might be due to the fact of CS is not only enhanced in regions with 
strong UV radiation, but its emission also increases towards the densest molecular cores due 
to its high critical d ensity. It is interesting to note that the similar FUV field derived by 
Kramer et al.l (120051 ) towards the nucleus of M 83 and NGC 5194 might explain the different 
abundances observed in our sample mostly as a consequence of the dense molecular fuel de- 
pletion and the molecular cloud structure in their nuclei, and not so much by the differences 
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in the photodissociation fields affecting the nuclear ISM. However , the CS variations in the 
Galactic center clouds are more moderate (a factor of ~ 2 — 3, iMartm et al.l |2008| ) than 
those of HNCO, and therefore it is likely to be averaged out over the central several hundred 
parsecs covered by the single dish beam. 



4.3. The evolution of nuclear starbursts 



The highly contrasted variation of HNCO abundances over the galaxy sample in this 
work allows us to trace the dense molecular gas available for fueling the nuclear starburst 
of these galaxies. This is equivalent to defining a chemical timeline in the evolution of the 
starburst event that could be read from right to left hand side in Fig. [31 It is important to 
note that the time scales for the chemical proce ssess due to shocks and photodissociation 
by UV radiation (10'* — 10^ yr. iBergin et al.lll998l ) are extreme ly short as c ompared with the 
typical time scale of the starburst phenomenon (10'' — 10^ yr. ICoziollll996l ). 



In this scenario, at the onset of the starburst, large amounts of dense molecular material 
are driven and compressed into the nuclear region of the galaxy unleashing the star formation. 
High abundances of HNCO are injected into the gas-phase within these clouds due to shocks, 
likely as a consequence of the clo ud- cloud collisions, as observed in the nuclear spiral arms 
of IC342 (IMeier fc Turner! |2005| ) . From our galaxy sample, M83 is shown to have the 
highest abundances of HNCO, and similar to those observed in the Galactic center GMCs. 
Subsequently, as soon as the first newly massive star clusters are formed, UV illumination 
begins to come on stage, photodissociating the material surrounding them. The HNCO 
abundance in these regions will be quickly dissociated where the CS molecule will be formed. 
A sort of average equilibrium will be reached between the PDR and the dense gas contribution 
as long as the infall of molecular fuel keeps feeding the starburst. This is the status in which 
most of the starburst galaxies in our sample are observed. The AGN dominated galaxies 
we have observed also show this average CS/HNCO abundances similar to the intermediate 
state starbursts, suggesting the presence of a significant UV radiation field. Once the galaxy 
gets closer to the end the short few 10 Myr starburst period, the fuel supply is cut off or 
stopped by superwinds. From our comparative study, this is apparently the point at which 
NGC253 stands, where photodissociation becomes more important and the fuel reservoir 
is being consumed by the last gasp of star formation. Finally, as we observe in M 82 and 
NGC5194, most of the molecular gas is pervaded by the UV radiation field with barely no 
dense UV shielded molec ular material left. The galaxy remains then waiting for a likely next 
burst of star formation ( Coziol "l996h. Indeed , chemical enrichment in M82 points out to 



recursive starburst periods (lOriglia et al.l 12004 ). 
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5. Conclusions 

From our study of the CS and HNCO abundances in a selected sample of extragalac- 
tic sources, we liave found differences of nearly two orders of magnitude in the HNCO/CS 
abundance ratios among the starburst galaxies. These differences are in agreement with the 
results from the observations of Galactic center clouds dominated by various heating mecha- 
nisms suggesting that the observed HNCO abundances could be related to the evolutionary 
stage of their nuclear starburst. Shocks affecting the massive dense molecular clouds in the 
early stages and strong UV fields in the evolved stages dominate the HNCO abundance in 
galaxies. We do not find significant differences in the molecular abundances of HNCO in 
galaxies with different type of nuclear activity (i.e. AGN or SB) which implies this species 
are not affected differently by these processes at the scales observed by single dish beam 
telescopes. 

As derived from the HNCO abundances we can point out prototypes of these nuclear 
starburst evolutionary stages. M83 and the nuclear spiral arms in IC342 would be repre- 
sentative of the early stage, while M82 and NGC5194 nuclei would be in the later stages of 
evolution. NGC253, claimed to be a prototype of early starburst, shows abundances below 
the average suggesting an important influence of photodissociation in the heating of its ISM, 
and therefore chemically closer to M82 than to M83. 

The ^^S/^'^S isotopic ratio, derived from the C'^^S/C'^^S abundance ratio, is found to be 
~ 10, and thus similar for all the observed sources. This ratio, found to be well below the 
value of ~ 24 in the Galactic disk, supports the idea of the ^"^S isotope being enriched due 
to massive star formation events in the nuclear region of galaxies. 

The CS emission is suggested to be a good tracer o f the molecular gas affected by 



photodissociation in the Galactic center (IMartin et al.ll2008l ). We do not observe a significant 



contrast in our sample. Though its abundance is enhanced in UV radiated regions, its 
emission is also favored in high density molecular clouds. Thus the overall emission over the 
central few hundred parsecs of galactic nuclei is averaged out, and no significant trends in 
the CS abundances is found in the observed galaxies. 
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Fig. 1. — Observed spectra for each source in the sample. Some of the transitions C^^O 1 — 
and C^^S 3 — 2 spectra have been scaled up for the sake of clarity in the profile comparison. 
Spectra with line intensities below 10 mK have been resampled to a velocity resolution of 
20kms^^, otherwise original velocity resolution is shown. C^^S (3 — 2) in NGC 1068 and 
C^'^S (2 — 1) in Arp220 are shown with a 30 km s~^ velocity resolution. Gaussian fits are 
shown for all detected transitions. Dotted fits indicate individual velocity components. For 
Arp220, where C^^O and HNCO transitions appear significantly blended, we show both fits 
in the same spectrum where the overall fit is shown with a continuous line and the individual 
transition with dotted lines. Temperature scale is Tub (mK). 
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Fig. 2. — C'^^S/C'^^S Abundance ratio for the detected galaxies as a function of the CS 
fractional abundance. Sulfur isotopic ratios from the literature for the Galactic Disk, and 
the galaxies NGC 4945 and NGC 253 are represented with continuous, dotted and dashed 
lines, respectively, and their corresponding error shown as a grey shade. Labeling of sources 
is the same as in Fig. [3l 
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Fig. 3. — Fractional abundances of CS {Top) and C^'^S {Bottom) vs the fractional abundance 
of HNCO with respect to H2. Upper limits to the abundances are represented by arrows. As 
in the text, offset observed position are labeled with an * following the galaxy name. The 
two velocity components fitted in Maffei2 are labeled as Maffei2a and Maffei2b. Dashed lines 
represent the limits of the regions used to discriminate the UV dominated (left side) from 
the shock dominated (r ight side) molecular gas, as derived from Galactic center observations 



see 



Martin et al.l 120081 . for details) 
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Table 1: Observational parameters 





SOURCE SAMPLE 




Source 


«J2000 


(^72000 


Dl/a (Mpc-pc/'O - 


NGC253 


00:47:33.3 


-25:17:23 


2.5 - 12 


Maffei2 


02:41:55.2 


59:36:11 


5.0 - 24 


NGC1068 


02:42:40.9 


-00:00:46 


16.7 - 81 


IC 342 ^ 


03:46:48.6 


68:05:46 


3.7-18 


M82 


09:55:51.9 


69:40:47 


3.5 - 17 


NGC4945 


13:05:27.2 


-49:28:05 


7.8 - 38 


NGC5194 (M51) 


13:29:52.7 


47:11:43 


9.7-47 


M83 


13:37:00.9 


-29:51:56 


3.7-18 


Arp 220 


15:34:57.1 


23:30:12 


73.0 - 354 


NGC6946 


20:34:52.4 


60:09:14 


5.5-27 


NGC 7469 


23:03:15.6 


08:52:26 


66 - 320 




MOLECULAR LINES 




Molecule 


Frequency (GHz) 


ft ("\ "5 

^boam ) 


f 

-Deflf 


Ci^O J = 1 - 


109.782 


22 


0.75 


Ci^O J = 2 - 1 


219.560 


11 


0.55 


CS J = 5 - 4 


244.935 


10 


0.49 


C^^S J = 2 - 1 g 


96.412 


25 


0.76 


J = 3 - 2 


144.617 


17 


0.69 


HNCO 5o,5 - 4o,4 


109.905 


22 


0.75 


HNCO 61,6 - 5i,5 


131.885 


19 


0.71 


HNCO 10o,io - 9o,9 


219.798 


11 


0.55 



Note. — The data for sources in italics were taken from previously published data. 
''Distances as derived bv lBaan et all (j2008l ) and linear scales given in parsecs per arcsecond at those distances. 
^Central and (0", -16") offset (NGC 1068* in the text) positions were observed. 
'^Central and (5", 15") offset (IC342* in the text) positions were observed. 

'^Central position was not observed, but the NE molecular lobe at the approximate offset (13", 7.5") position, 
referred to as M 82* in the text. 

^All data are from the IRAM 30 m telescope, but for those of NGC 4945 from SEST |Wang et al.ll2004h . In 
this case, the beam is about twice the given size. 
^Beam efficiency of the IRAM 30 m at the observed frequencies 
^This transition is only observed towards Arp 220. 
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Table 2: Parameters derived from the Gaussian profiles fitted to the observations. 



Transition 




"LSR 


^"1/2 




rms 




(Kkms^l) 


(kms^l) 


(kms^-') 


(mK) 


(mK) 






Maffei 2 








C^^O 1-0 


1.2 ± 0.2 


-102 ± 4 


63 ± 7 


18.2 


1.9 




3.7 ± 0.3 


4 ± 4 


120 ± 10 


29.0 




C^**0 2-1 


2.0 ± 0.2 


-102.2 ± 1.4 


66 ± 3 


28.0 


1.0 




8.0 ± 0.2 


-14 ± 2 


151 ± 4 


49.8 




CS 5 - 4 


1.1 ± 0.2 


-93 ± 4 


63 ± 9 


16.3 


1.4 




1.7 ± 0.3 


2 ± 10 


110 ± 20 


14.3 




C34g 3 _ 2 


0.41 ± 0.08 


-100 


65 


6.0 


1.4 




0.37 ± 0.12 





110 


3.2 




HNCO 5—4 


0.38 ± 0.10 


-100 


62 


5.8 


2.1 




1.42 ± 0.13 


10 ± 3 


66 ± 7 


20.1 




HNCO 6-5 


< 0.3 




65 


< 5.0 


2.5 




1.67 ± 0.17 


4 ± 3 


66 ± 8 


23.9 




HNCO 10 — 9 


0.13 ± 0.06 


-100 


65 


1.8 


1.0 




0.39 ± 0.06 


5 ± 6 


62 ± 10 


5.9 






NGC 1068 








C O 1-0 


4.45 ± 0.12 


1102 ± 4 


264 ± 8 


15.8 


0.8 


C-'^^O 2—1 


2.0 ± 0.2 


1067 ± 12 


210 ± 30 


9.0 


1.9 


CS 5 — 4 


3.3 ± 0.3 


1125 ± 9 


180 ± 20 


17.6 


3.0 


C'^'^S 3 — 2 


0.53 ± 0.12 


1070 ± 30 


230 ± 60 


2.2 


1.0 


HNCO 5 — 4 


0.86 ± 0.09 


1086 ±11 


210 ± 20 


3.9 


0.8 


HNCO 6 — 5 


0.56 ± 0.07 


1114 ± 13 


210 ± 30 


2.5 


t).6 


HNCO 10 — 9 


< 0.4 




210 


< 2.1 


1.9 




NGC 1068 (O", -16") 








IS 

C O 1-0 


3.39 ± 0.14 


1130 ± 4 


213 ± 10 


15.0 


1.1 


C^^S 3-2 


< 0.2 




170 


< 1.4 


1.2 


HNCO 5 — 4 


0.41 ± 0.11 


1160 ± 20 


170 ± 40 


2.3 


1.1 


HNCO 10 — 9 


< 0.2 




170 


< 1.0 


0.8 




IC 342 (()", O") 








IS 

C O 1-0 


4.80 ± 0.11 


25.2 ± 0.6 


50.1 ± 1.2 


89.9 


1.7 


C^^O 2—1 


8.10 ± 0.10 


29.9 ± 0.3 


46.3 ± 0.6 


164.1 


1.1 


CS 5 — 4 


2.15 ± 0.11 


37.7 ± 1.2 


50 ± 3 


40.3 


1.8 


C^'^S 3 — 2 


0.49 ± 0.12 


26 ± 7 


54 ± 12 


8.4 


1.5 


HNCO 5 — 4 


1.81 ± 0.10 


30.5 ± 1.4 


52 ± 3 


33.0 


1.7 


HNCO 6 — 5 


1.87 ± 0.11 


35.0 ± 1.6 


53 ± 3 


33.0 


1.9 


HNCO 10 — 9 


0.31 ± 0.06 


24 ± 3 


34 ± 8 


8.5 


1.1 




IC342 (5", 15") 








s 

C O 1-0 


3.45 ± 0.10 


43.9 ± 0.6 


41.1 ± 1.4 


78.8 


2.0 


CS 5 — 4 


0.42 ± 0.10 


49 


36 


11.0 


2.6 


C^^S 3-2 


0.32 ± 0.04 


49 ± 2 


36 ± 6 


8.3 


1.0 


HNCO 5 — 4 


1.39 ± 0.10 


44.9 ± 1.2 


36 ± 3 


36.2 


2.0 


HNCO 6 — 5 


0.81 ± 0.10 


41 ± 2 
M 82 


34 ± 5 


22.7 


2.0 


s 

C O 1-0 


1.55 ± 0.07 


235 ± 2 


87 ± 5 


16.8 


0.7 




4.41 ± 0.06 


310 


80 


51.9 




C *0 2—1 


22.79 ± 0.15 


309.5 ± 0.3 


81.9 ± 0.6 


261.5 


2.0 


CS 5 - 4 


6.4 ± 0.2 


313.0 ± 1.0 


73 ± 2 


82.6 


2.5 


C34g 3 — 2 


0.58 ± 0.08 


302 ± 6 


91 ± 13 


6.0 


1.0 


HNCO 5 — 4 


< 0.10 




80 


< 1.3 


0.7 


HNCO 6 — 5 


< 0.17 




80 


< 2.3 


1.3 


HNCO 10 — 9 


< 0.3 




80 


< 3.7 


2.0 




NGC 5194 








IS 

C O 1-0 


0.80 ± 0.08 


418.2 ± 1.4 


39 ± 3 


19.1 


0.6 




0.86 ± 0.09 


490 ± 5 


80 ± 8 


10.1 




C S 3 — 2 


< 0.2 




100 


< 2.2 


1.4 


HNCO 5—4 


< 0.11 


M 83 


100 


< 1.0 


0.6 


s 

C O 1-0 


2.19 ± 0.14 


499 ± 2 


66 ± 5 


30.9 


1.9 


CS 5 — 4 


2.3 ± 0.3 


504 ± 7 


82.4 ± 1.1 


26.5 


0.5 


C^'^S 3 — 2 


0.40 ± 0.13 


501 ± 6 


34 ± 13 


11.1 


2.7 




0.70 ±0.15 


510 ± 10 


90 ± 19 


7 7 


1 9 


HNCO 6-5 


0.38 ± 0.07 


515 ±7 


67 ± 12 


5.4 


1.2 


C^^O 1-0 




Arp 220 








2.7 ± 0.3 


5445 


450 ± 40 


5.7 


0.9 


CS 5 - 4 


11.3 ± 0.2 


5446 ± 5 


481 ± 11 


22.0 


2.7 


C^^'S 2 - 1 


0.9 ± 0.2 


5460 ± 30 


240 ± 50 


3.4 


1.4 


HNCO 5-4 


1.1 ± 0.2 


5445 


450 


2.3 


0.9 




NGC 6946 








C^^O 1-0 


1.36 ± 0.10 


-5.0 ± 1.4 


54 ± 3 


23.8 


0.7 




3.41 ± 0.11 


84.7 ± 1.7 


103 ± 4 


31.2 




C^*S 3-2 


0.37 ± 0.07 


71 ± 11 


90 ± 20 


3.8 


0.7 


HNCO 5-4 


0.55 ± 0.10 


62 ± 8 


87 ± 19 


5.9 


0.7 


HNCO 6-5 


0.24 ± 0.06 


54 ± 7 


53 ± 14 


4.1 


1.0 




NGC 7469 








C^^O 1-0 


0.61 ± 0.05 


4910 ± 10 


250 ± 20 


2.3 


0.4 


0^**3 3-2 


< 0.13 




150 


< 0.9 


0.6 


HNCO 5-4 


0.16 ± 0.04 


4890 ± 20 


150 ± 40 


1.0 


0.4 
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Table 3: Derived column densities and rotational temperatures for each galaxy. 



Source 


Ci8 





CS 




HNCO 


( 


xlO^^cm-^) 


(K) 


(xlO^^cm-^) 


(xl0i2cm-2) 


(xlO^^cm-^) 


(K) 


NGC253 ^ 


100 ±20 


5.6 ±0.4 


21 ±4 


19 ±4 


22 ±9 


18 ±4 


Maffei 2 ^ 


5.7± 1.5 


5.6 ±0.6 


2.4 ±0.7 


5.8 ±1.1 


2.5 ± 1.2 


12.2 ± 1.9 




17 ±2 


6.6 ±0.3 


3.5 ±0.9 


5.2 ± 1.0 


10 ±2 


11.5 ±0.6 


NGC1068 


24 ±2 


3.3 ±0.1 


7.0 ±1.5 


7.3 ± 1.2 


10 ±5 


5.8 ±0.8 


NGC1068* 


18 ±4 






< 3.4 


2.9 ± 1.0 




IC 342 


22 ± 1 


5.7±0.1 


4.6 ±0.1 


6.9 ± 1.3 


12 ± 1 


10.2 ±0.4 


IC 342* 


18 ±4 




0.9 ±0.2 


4.4 ± 1.0 


19 ± 1 


5.3 ±0.6 


M82* 


28 ± 1 


14.5 ±0.3 


14 ±4 


8.1 ± 1.9 


< 0.7 




NGC4945'= 


160 ± 10 


6.5 ±0.2 


58 ± 12 


64± 13 


50± 16 




NGC5194 


8.6 ± 1.8 




< 0.7^^ 


< 3.0 


< 0.7 




M83 


11 ±2 




4.9 ± 1.2 


5.6 ± 1.1 


10 ±5 


5.0 ± 1.1 


Arp 220 


14 ±3 




24 ±5 


26 ±5 


7±3 




NGC6946 


25 ±5 




6.0 ±0.8^ 


5.2 ± 1.0 


12 ±6 


4.2 ±0.9 


NGC 7469 


3.2 ±0.7 






< 1.9 


1.1 ±0.4 





* Indicates offset positions. 

" Derived using data from [Harrison et all |l999l) . iMartfn et all |2005l ): iMartfn et all |2006bl ). 
Abundances derived for each of the two velocity components labeled as Maffei 2a and Maffei 2b in Fig. [Sj 
Derived using data from lWang et al.l (j2004l ). 
Derived using data from Mauersberger et al. (198^. 



